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The splitting of matter-waves into a superposition of spatially separated states is a funda- 
mental tool for studying the basic tenets of quantum mechanics and other theories, as well 
as a building block for numerous technological applications. We report the realization of a 
matter-wave beam splitter based on magnetic field gradients on an atom chip, which can 
be used for freely propagating or trapped atoms in a Bose-Einstein condensate or a thermal 
state. It has a wide dynamic range of momentum transfer and operation time. Differential 
velocities exceeding 0.5 m/s can be achieved in a few micro-seconds. The beam splitter may 
enable a wide range of applications, such as, fundamental studies of many-body entangle- 
ment and dephasing processes, probing classical and quantum properties of nearby solids, 
and metrology of rotation, acceleration and gravity on a chip scale. 



I. INTRODUCTION 

For the last two decades, matter- waves, in the form of ultra cold atoms [1] [2], have been the 
source of numerous new experimental insights. One of the basic tools for such studies is matter- 
wave interferometry which enables high precision measurements for fundamental physics [SHE] as 
well as for technological applications [9HT3]. 

In one type of matter-wave interferometry, atoms are split into a superposition of momentum 
states by interaction with laser light, where the signal is typically derived from the population 
of internal states |14H18| . High precision measurements, e.g. of acceleration and rotation, are 
enabled by an accurate momentum transfer from laser photons to atoms in multiples of the photon 
momentum hk, where k = 2ir/\ is the wave- vector of the photons. 

In another type of interference experiments, a Bose-Einstein condensate (BEC) is split into two 
parts by a potential formed by optical \19\ I20j. radio-frequency (RF) or micro- wave \21\ I22j fields. 
The signal is typically derived from spatial interference patterns. These experiments provide new 
experimental insight in topics such as coherence and entanglement in a many-body system |23H25| . 
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Potential applications of spatial interferometry with a BEC on an atom chip [26- 28j , which is the 
setting used here, include the study of properties of many body systems, as well as the sensitive 
probing of classical or quantum properties of solid state nano-scale devices and surface physics 
(e.g. [29H33]). The latter is expected to enhance considerably the power of non-interferometric 
measurements with ultracold atoms on a chip, which have already contributed, for example, to the 
study of long-range order of current fluctuations in thin films |34| , the Casimir- Polder force |35H38| 
and Johnson noise from a surface [39j. 

In addition, coherent splitting of matter-waves on a chip is a crucial step towards the develop- 
ment of miniature rotation, acceleration and gravitational sensors based on guided matter-waves 
|4CH I41j. Depending on the progress made in developing stable electronics, a current-carrying atom 
chip may be used for metrology with either guided or freely propagating atoms. 

Here we propose an atom chip field-gradient-beam-splitter (FGBS), which may provide an 
efficient and flexible tool for interferometry in the micrometer scale as well as a macroscopic scale. 
The operation of the FGBS is simple and in the present realization requires only the manipulation 
of static magnetic and RF fields. The FGBS enables a large momentum transfer in a short time, 
and as the coherence time of split matter-waves is limited by atom-atom interaction and coupling 
to the environment, it may enable sensitive interferometry, for example by way of large density 
and large area, within the available coherence time. We present a specific realization of the FGBS 
using magnetically sensitive Zeeman levels and demonstrate the formation of spatial interference 
fringes as an interferometric signal for both freely propagating and trapped atoms. We also propose 
other realizations using first-order magnetically insensitive states ("clock states") and internal state 
population as a signal. 



II. GENERAL SCHEME 



The general scheme of the FGBS and its output is demonstrated in Fig. [T} It is based on 
a Ramsey- like sequence of two tt/2 rotations with a field gradient applied during the interval 
between them. If the gradient is applied perpendicular to the motion of the atoms, the differential 
momentum will be in the transverse direction and an area enclosing interferometry would ensue. 
In our experimental demonstration the gradient is applied along the direction of motion to form 
a one dimensional interferometer. Considering two-level atoms with internal states |1) and |2) 
and starting with atoms in state |2), the first tt/2 rotation transfers them into the superposition 
state -4=(|1) + |2)). We then apply a field gradient which constitutes a state-selective force Fj(j = 
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FIG. 1. Absorption images of the atoms and the corresponding experimental sequence: (a) in the trap 
before the release, (b) after weak splitting of less than hk in 5/xs interaction time and 14 ms time-of-fhght 
(TOF), (c) after a strong splitting of above 40 hk in 1ms interaction time and 2 ms TOF, and (d) the four 
clouds separated by another strong gradient (Stern-Gerlach) after the FGBS. Images (b) and (c) show the 
high dynamic range this method can reach without any complicated sequence. Image (d) shows that each 
cloud in (b) and (c) is constructed from atoms in the |1) and |2) states (see text for details). 

1,2) = — VVj- for an interaction time T, which is typically shorter than the time it takes the atoms 
to move in the force field. The potential V may be due to any kind of field such as electro- magnetic, 
electric or magnetic, and the states |1) and |2) may be any states between which one can introduce 
a coupling and for which there exists a state selective force. The state of the atoms after time 
T is given by |^(x,T)) = X (|l)e iFrxT/ri + |2)e iFVxT/?l ) V'o(x), where ^o(x) is the spatial wave- 
function at t = 0. Each level gains a phase </>(x, t) = —AEjt/h = Fj-xT/h, which is equivalent to a 
momentum transfer pj = FjT. The second 7r/2 rotation transfers the atoms into the superposition 
state 

|V) = \ (|l, P i) + |1, P2 ) - |2,pi) + |2, P2 » (1) 

such that each of the internal states \j) is in a superposition of two momentum states \pj) = 
e* Pj " x /^o( x ), and vice versa. Eq. (1) is also valid if the atomic motion during the interaction 
time is taken into account or when the force is not homogeneous in space and time, given that the 
momentum states \pj) are more generally defined as the states that evolve from ipo(x) under the 
influence of the respective potentials Vj(x,t). 

This completes the operation of spatial splitting. One may then use another field gradient pulse 
to separate between the |1) and |2) states [Fig. [T|(d)] , and use these states to realize two parallel 
interferometers for noise rejection. One can also use the entangled momentum and internal state 
as an interferometer of two clocks (e.g. [12] )• If one wishes to stay with just one of the two internal 
states, one may typically find a dedicated transition to discard the redundant state. 
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III. EXPERIMENT AND RESULTS 

To realize the FGBS, we utilize Zeeman sub- levels. The idea of splitting different spin states 
in a magnetic gradient is as old as the Stern-Gerlach experiments of the 1920s (e.g. [12]), but to 
the best of our knowledge was never designed with the necessary detail and to this day remained 
unrealized. Furthermore, unlike such schemes, the FGBS gives rise at its output to two wavepackets 
with identical spin states, an advantageous feature in noisy environments. 

A. Splitting free falling atoms 

We start with a BEC of ~ 10 4 87 Rb atoms in state \F, m F ) = |2,2) = |2) and use an RF field 
to perform transitions to state |1) = |2, 1). We utilize the same set-up described in [Hj. The trap 
position is 100 /iin from the chip surface, and the radial trapping frequency of 1 2) is 002 ~ 2-7T x 100 Hz 
(for state |1), u)\ = U2/V2). In order to have the |1) and |2) states form a pure two-level system, 
we apply a strong homogeneous magnetic field (A-EJ12 ~ h X 25 MHz) and push the transition to 
|2,0) out of resonance, due to the nonlinear Zeeman effect (~ 250kHz) [Hj. Next, we release the 
BEC and apply two n/2 RF pulses [45] with a Rabi frequency of = 20 — 25 kHz and with a 
magnetic gradient pulse of length T in between, thus forming a Ramsey-like interferometer. 

This scheme for freely falling atoms has been applied for atoms above and below the BEC 
transition temperature. Here we present our results for a BEC. Fig. [ija) presents the cloud 
position in the trap, before it was released. In Fig. [l][b,c) we present two splitting events which 
exhibit the large dynamic range of this method, one with a differential momentum of less than hk 
(T = 5/is), and the other with above 40 hk (T = 1ms). Here hk is a reference momentum of a 
photon with 1 jim wavelength. The only parameter which is changed is the interaction time T, in 
which a current of 2 — 3 A in a 200 x 2 fim 2 gold wire on the chip surface gives rise to the gradient. 
To verify the internal state of the atoms [Eq. ([!])] we separate between the |1) and |2) states by a 
long pulse of magnetic field gradient (Stern-Gerlach), as shown in Fig. [T][d). 

Fig. [2] shows the measured differential momentum transfer as a function of the interaction time 
T. An interaction time as short as 100 /xs is required to transfer a relative velocity of 50mm/s 
(equivalent to 10 hk). Continuing the general theoretical description given above, the operation 
of the magnetic realization of the FGBS may be explained by simple kinematics. During the 
interaction time T, a differential acceleration between the wavepackets is induced, such that after 
the FGBS, each internal state is a superposition of two wavepackets which were accelerated as a 
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FIG. 2. Characterization of the FGBS for free falling atoms: The differential velocity between the two 
observed clouds as a function of the interaction time T. The error bars are calculated from the variation of 
a few data sets. The theoretical prediction (shaded area) is based on a measured current in a resistor that 
mimics the chip wire [inset (a)], taking into account errors of ±2.5 /mi (one pixel) in the initial cloud position, 
±0.2 O in the wire resistance and a 1 us delay of the measured rise-time (due to the resistor's inductance). 
The linear dependence on T is to be expected from the solution of Eq. The current "overshoot" at 
short times is responsible for the larger acceleration for small T. Inset (b) : Numerical integration of Eq. ^ 
over T, using the experimental wire configuration, and using a constant current of 3 A. We observe the 
limit of our method due to the large distance which develops between the atoms and the wire. Inset (c): 
Utilizing improved experimental parameters we find that Ap of above 100 Hk is possible in less than 10 /is. 
The parameters are: z = 10 |xm, / = 2 A and the wire dimensions are 10 x 2/xm 2 (10 7 A/cm is more than 
sustainable for such short pulses). 



|1) or |2) state. The momentum kick for a wavepacket of a certain tuf state at a distance z below 
the chip wire is 



dp 
~dt 



2ttz 2 



(2) 



where fio and \xb are the magnetic permeability of free space and the Bohr magneton, gp is the 
Lande factor for the hyperfine state F, and I is the current. The equation does not present the 
nonlinear term in B and a geometric term 1/[1 + (W/2z) 2 ], accounting for the finite width W of 
the wire. These terms have not been neglected in our full simulation of the FGBS (in Fig. [2J which 
predicts the experimental results and shows that momentum transfers of over 100 hk are feasible 
in a few micro-seconds [see inset (c), Fig. |2j. The linear relation in Fig. [2] is to be expected for the 
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FIG. 3. Simulation of the free fall experiment: For (a) T = 5 /is and (b) T — 10 /xs, we present the atomic 
density of the |2) state (GP simulation), just after the splitting. These "near field" fringes are Fourier 
transformed into two wavepackets in the "far field". This point of view is an alternative to the simple 
kinematic explanation provided by Eq. ([2]). 

short interaction times during which the atoms move only slightly and the acceleration in Eq. ^ 
is fairly constant (the kink in Fig. [2] is due to changing currents, see caption). 

In Fig. [3] we present a Gross-Pitaevskii (GP) simulation of the atomic density for state |2) just 
after the second ir/2 pulse ("near field"). The density shows a sinusoidal pattern. This pattern 
may be viewed as spatial Ramsey fringes: an atom at position x is subjected to a Ramsey sequence 
in which the phase accumulated between the two ir/2 pulses is 1/h AV(x, t) dt where AV is 
the potential difference and T is the interaction time. As the potential varies linearly along the 
z axis, Ay = AFz, and at short times the atom stays almost stationary, the Ramsey phase 
and consequently the population of the two states |1) and |2) is modulated with a wavelength 
A = h/ AFT. The formation of sinusoidally modulated internal state populations is equivalent 
to varying Rabi frequencies along the ensemble when the driving field has an intensity gradient 
|46H48j . when AF is replaced by the gradient of the Rabi frequency (up to K). As a perfect sine 
function Fourier transforms into two identical and counter propagating k components, this "near 
field" density distribution transforms after a sufficient time-of-flight into a pair of wavepackets with 
momentum difference Ap = h/\. We note that this is the inverse of a recombination process of 
two wavepackets. 



B. Splitting trapped atoms 

Next, we show that our scheme could also be implemented on a trapped BEC, thereby opening 
the road for guided "multi-pass" interferometry, such as the Sagnac we have previously proposed 
|4U| . We repeat the above scheme and performing all the sequence on an atom cloud trapped in 
an Ioffe-Pritchard (IP) magnetic trap. We start with ~ 10 4 87 Rb atoms in a trap 250 /im from 
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FIG. 4. Splitting a trapped BEC: (a) An illustration of the two traps below the atom chip showing the 
vertical separation between the two traps, and the difference in size due to the different confinement, (b-d) 
Trapping potential of the two states (the energy difference was minimized for visibility) in the direction 
of gravity during the splitting process. The four parts of the wave-function and their momentum (size of 
arrow) are also shown, (b) The splitting just after the first 7r/2 pulse, where the position of both clouds is 
the trap minimum of the |2) state, (c) The clouds after the second ir/2 pulse where only the |1) part of the 
super-position was accelerated before the pulse. There is still spatial overlap (the distance in exaggerated 
in the image for clarity), (d) The four parts after some oscillation time in the trap, showing that one part 
of the 1 2) state didn't move, the other part of the |2) was slowed by the trapping potential almost to a halt, 
while the two parts of |1) were accelerated. 

the chip, with a radial trapping frequency of U2 ~ 2ir x 100 Hz. Due to the high magnetic field 
at the trap bottom, AE12 ~ h X 18 MHz and the nonlinear Zeeman shift is about 100 kHz. The 
sequence is similar to the free falling experiment but with the release taking place after the second 
7r/2 pulse. Contrary to the free falling case, no external gradient is required as the trap gradient 
naturally accelerates the atoms. Furthermore, the gradient exists also during the it/ 2 pulses, but 
due to the wide spectrum of the short pulse it is still quite effective in creating the superposition 
(a 7r pulse with tin = 5 — 10 kHz transfers up to 90% of the atoms, compared with more than 95% 
in free fall). 

In Fig. Qa) we present the two potentials by their equipotential surfaces, as they are situated 
below the atom chip. In Fig. |4fb-d) we utilize a ID energy versus position (z) plot to describe 
the evolution of the system during and after the FGBS sequence. Due to gravity, the centers 
of the combined magnetic and gravitational trapping potentials for the two levels are shifted by 
Az = g/u)2 2 - It follows that when atoms initially at the level |2) are transferred by the first ir/2 
pulse into level |1), they experience an acceleration 

a(T) = bj\ 2 z{T) = ufAz cos(u}\T) = ^cos(wiT). (3) 

For interaction times T <C tt/uj\ these atoms move only sightly along the potential gradient such 
that the acceleration is constant and equals g/2, and as in the free fall scheme the momentum 
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FIG. 5. Characterization of the FGBS for trapped atoms: The differential velocity of the two wavepackets 
as a function of the interaction time T. Av, and its error, were taken from a linear fit of six measurements 
for each T, at two different Rabi frequencies and three TOF. As expected from Eq. Av grows linearly 
with time. The solid line is a theoretical curve v(T) = | sin(wiT)/aji + v r , where v r = 0.58 mm/s is an 
additional velocity due to atom-atom repulsive interaction (no fitting parameters). The first term follows 
from Eq. Q , while the collisional constant v r is obtained from a full GP simulation. The linearity of the 
graph is due to the small interaction times such that siti(ujT)/uj w T. 

splitting grows almost linearly with T. This almost linear dependence is shown in Fig. [5j Here, 
the atom-atom collisional repulsion is responsible for an additional contribution to the velocity, as 
shown in the experimental data and confirmed by our numerical GP solution. 



C. Recombining the two wavepackets 



In Fig. |6^a,b) we present typical images of interference fringes as observed in our free-fall 
experiments, where we use an additional gradient pulse at a specific time At after the first one. 
At this time, the two wavepackets are already separated by a distance d = ApAt/m, Ap being the 
momentum transfer difference due to the FGBS. The second gradient pulse gives a stronger kick 
to the wavepacket with smaller momentum, which is at this time closer to the chip wire relative 
to the wavepacket with larger momentum. The duration of the second momentum kick is tuned 
such that after this kick the two wavepackets have the same momentum (with spatial separation 
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FIG. 6. Interference fringes: (a) Large (T = 5 /is) and (b) small (T = 10 /is) periodicity fringes, created 
by recombining the two outputs of the FGBS. Cuts of the optical density (OD) are also provided, with 
fits (see Methods for fit description). The two |2) wavepackets were recombined by an additional magnetic 
gradient and imaged after 14 ms of TOF. The different T gives rise to a different At>, and hence a different 
distance before recombination, leading to different fringe spacings, 33 /im and 16 /im, respectively, (c) Fringes 
observed after the FGBS is applied to a trapped BEC (T = 700 /is, TOF = 10 ms), when recombining the two 
1 2) wavepackets by allowing them to oscillate 2 ms in the trap before releasing them. The 8 /xm periodicity 
fits well with the well known formula for fringe separation ht/md. The low visibility is attributed mostly to 
our limited optical resolution (ss 5 /jm) and our less-than-perfect ir/2 pulses. 

d). They expand and overlap after a TOF to give rise to an interference pattern, which is fitted to 



where zq and a are the center and length of the joined wavepacket, A = ht/md and (ft are the 
wavelength and phase of the fringes, and v is the visibility. 

Fig. [6^c) shows an interference pattern generated with a BEC split in the trap. Here we have 
allowed the atoms to oscillate in the trap for a period of about 2 ms, which is approximately a 
quarter of the trap harmonic period, so that the relative velocity between the two wavepackets is 
almost zero, see Fig. [4](ci) . The trap is then released and the two wavepackets, positioned at z ~ 
and at z ~ Ap/muj = d, expand, overlap, and form multiple interference fringes. The fit in Fig. 
ptc) used a form similar to Eq. Q but with a Thomas-Fermi envelope 







and the data was actually fitted to the form 
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which takes into account a reduction due to saturation. 

IV. PHASE AND MOMENTUM STABILITY 

Beyond the presentation and characterization of the technique and kinematics of the FGBS, it 
is important to analyze the prospects and bounds on the stability of the phase and momentum 
differences, A(j) and Ap, between the two outports of the FGBS. These stabilities determine the 
quality of interferometry that the FGBS can support. For a general FGBS, 

A<t>=(y 2 -V 1 )T/himd Ap=(F 2 -F 1 )T, (7) 

and the stability is derived from the specific implementation of the FGBS. For magnetic fields 
generated from thin chip wires and considering only the linear Zeeman effect (with our |1) and |2) 
states), the FGBS imparts 

IT IT 

A6 oc — and Ap oc (8) 

z z z 

and consequently 5(p/A<f> and dp/Ap, the relative instabilities at the output of the FGBS, are 
proportional to SI /I, 5T/T and bzjz. Our main source of instability, which led to 5(f) ~ 27r, was 
the latter term, as the trapping position was not determined by the gradient inducing wire but 
rather by a distant wire. 

The position instability can be made negligible by a trapping procedure based on a single current 
such that the trap position depends only on the geometry and is independent of the current. For 
example, let us use, for both the trapping and the momentum kick, three parallel chip wires 
(distance D) connected in series, where the two side wires replace the external bias field. The trap 



position z = D is independent of I. As B 



2tt 



1 2z 



z z 2 +d 2 



and 5B = f B + §f 5z, the small Sz 



will suppress the second term in 5B. In addition, the created quadrupole will suppress the first 
term as well, as it has the same gradient as a single wire but with a much weaker magnetic field. 
Consequently, 5(ft will be made to depend only on 5T as oc BT. 

In order to estimate the bounds on the 5p/Ap stability, one should consider the available 
technology. Let us assume a 10 [is pulse. Then, for a 2 A current (containing ~ 10 14 electrons), 
the shot noise leads to 51 /I ~ 10~ 7 . Power sources with sub-shot noise are being developed, e.g. 
at JPL |49j . and may enable an even better stability (one obviously needs to also account for 
flicker and Johnson- Nyquist noise). Ultra stable capacitor driven current pulses may also improve 
performance (at mK temperature stability, available capacitors reach 5C/C = 10 -9 ). For pico- 
second switching electronics, one similarly finds 5T/T ~ 10~ 7 . We note that for longer pulses, 
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the ST/T stability will improve linearly with time (e.g. ~ 1CT 9 for a one mili-second pulse), and 
thus with momentum transfer (as long as SI /I is not the limiting factor). Overall, a momentum 
stability of 10~ 7 and beyond, should be considered feasible. 

In the trapped BEC experiment, the main source of instability is different. The g/2 gradient 
is fixed and does not fluctuate, and because of the long time between the pulses, 5T/T can also 
be neglected. The main limitation is the shot-to-shot fluctuations of the magnetic field at the trap 
bottom, which changes the energy between the two levels. 

V. CONCLUSIONS AND OUTLOOK 

To conclude, we have presented a new beam-splitting scheme. The FGBS can make use of 
different types of potentials and internal states. The scheme utilizes a Ramsey-like process which 
introduces a differential acceleration between two wavepackets by a state selective gradient. In the 
specific realization presented here, the FGBS makes use of magnetic gradients and Zeeman sub- 
levels. The FGBS is shown to enable a very large splitting momentum in very short times while 
also allowing for a wide dynamic range. The new beam-splitter may be used for freely propagating 
atoms as well as a trapped BEC and may be utilized for future interferometry experiments both 
for fundamental studies and technological applications. 

This work focused on the spatial splitting of an atomic wavepacket and we have briefly presented 
an interferometric signal based on spatial fringes. As an outlook, and in order to further emphasize 
the versatility of the FGBS, let us note that while the present work made use of pulses in the time 
domain, a continuous beam apparatus may make use of pulses in the spatial domain. In addition, 
we consider the possibility of the FGBS to also give rise to a signal based on the population of 
internal states. Furthermore, we also show that different internal states may be used by the FGBS. 

A. Interference with population fringes 

Typical interferometers probe the phase between two propagation paths by employing a second 
beam splitter to translate this phase into populations of momentum states or internal atomic states. 
Such an interferometer based on the FGBS would include an appropriate "mirror" that will bring 
the two wavepackets back to the same point in space, each with a different momentum state. At 
this point, we apply another FGBS sequence that would give an appropriate differential momentum 
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kick. The quantum state of an atom just before the operation of the second FGBS is 

A=[\2,P+ Ap/2) + e**\2,p- Ap/2)), (9) 

where p = \{p\ + P2) is the average momentum just before the second FGBS, Ap is the differential 
momentum due to the first FGBS, and <f> is the phase difference between the two paths which was 
accumulated during the propagation. Based on Eq. ([!]), the state after the second FGBS, which 
again introduces a differential momentum Ap, would give 

^[cob(0/2)|1> +sin(^/2)|2)]|p / ) + (10) 
^=[(|1) + \2))\p f + Ap) + e **(|l) - |2»|p, - Ap)], 

where pf is the average final momentum after the interferometer. The phase <j> can then be 
determined by a selective detection of the population of the internal state of the atoms with 
the momentum state \pf). The drawback of this interferometry scheme is that half of the atoms 
are discarded in each beam-splitter and only 1/4 of the atoms contribute to the interferometric 
signal. 



B. Splitting magnetic insensitive states 



Another interferometric scheme which is based on the same idea of the FGBS involves the two 
magnetically insensitive sublevels |1,0) and |2,0) of the two hyperfine states, which are analogous 
to those used in present-day precision interferometers. A micro- wave n/2 pulse may be used to 
create an equally populated superposition of the two states and then they may be split into two 
momentum states using a magnetic gradient at a high magnetic field. The nonlinear Zeeman shift 
of the transition energy between the states is AE ~ aB 2 , where a = 2irh x 575Hz/G 2 . At a 
distance of 10 ^m from a wire carrying 2 A of current, the atoms are exposed to a magnetic field 
of B = 400 G and a magnetic gradient of d z B = 40kG/mm. It then follows that F = d z AE(z) = 
2nh x 575 x 2Bd z B = 1.22 x 10 -20 N and consequently atoms receive a differential velocity of 
84.7 mm/s for a 1 /js pulse, equivalent to about 18 hk. As the two states are relatively magnetically 
insensitive, a second n/2 pulse would not be needed and the two output beams of this beam splitter 
could be used for interferometry in a completely analogous way to existing interferometers based 
on light beam splitters, by, for example, employing a second tt/2 pulse when the wavepackets are 
recombined, and measuring the populations of the internal states. 
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